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Analytical and numerical calculations for a reduced Fermi-Pasta-Ulam chain demonstrate that en- 
ergy localization does not require more than one conserved quantity. Clear evidence for the existence 
of a sharp delocalization-localization transition at a critical amplitude A c is given. Approaching A c 
from above and below, diverging time scales occur. Above A c , the energy packet converges towards 
a discrete breather. Nevertheless, ballistic energy transportation is present, demonstrating that its 
existence does not necessarily imply derealization. 

PACS numbers: 05.45.-a, 05.60Cd, 63.20.Pw 



One of the classical investigations of relaxation dynam- 
ics of macroscopic systems is to determine the time evo- 
lution of a perturbed equilibrium state. If this initial 
state converges to equilibrium the system is called mix- 
ing, implying ergodicity, and it is nonmixing, otherwise. 
An important question is: Does there exist a sharp er- 
godicity breaking transition under variation of a physical 
control parameter like temperature or strength of pertur- 
bation? 

Within a mode coupling theory for supercooled liq- 
uids such a dynamical glass transition has been found, 
but its sharpness seems to result from the mode coupling 
approximations (for reviews see [H, 0| ) • It is n °t °ur pur- 
pose to contribute to the theory of glass transition, but to 
study the influence of anharmonicity on the relaxational 
behavior at zero-temperature. In that case the generic 
lowest energy state of a particle system is a crystal. One 
may ask a similar question as above: Does an initially lo- 
calized energy excitation spread over the complete crys- 
tal, or not? In case of small excitation amplitudes, one 
can use the harmonic approximation. Then the time 
evolution of an initial configuration can be determined, 
exactly For one-dimensional harmonic lattices the 
results are particularly simple 0, Q ■ Independent of the 
strength and size of the excitation it always spreads over 
the full system, and energy transportation is ballistic, 
provided that there is no disorder. That infinite har- 
monic crystals are ergodic in general, has been proven 
rigorously Q. If the excitation amplitude increases, an- 
harmonicity gets important. 

Let us neglect any disorder, but taking anharmonic 
interactions into account. Discreteness of the lattice 
combined with anharmonicity allows for the existence 
of localized periodic vibrations, called discrete breathers 
(DB). For reviews see Ref. [Tfl. Their existence sug- 
gests that under certain conditions a localized excita- 
tion may converge to a DB, whereby suppressing com- 
plete energy spreading. Indeed, numerical solutions of 
the discrete nonlinear Schrodingcr equation (DNLS) [8] 
and references wherein, the Klein Gordon chain (KG) 
@ and the /3-Fermi-Pasta-Ulam chain (FPU) pa, GJ 
demonstrate generation of DB and their role for slow en- 
ergy relaxation. Particularly, the numerical results in 
Refs. p| |Tl[ give evidence that DB generation from an 



initially localized excitation requires an excitation ampli- 
tude which is large enough. This has been supported by 
analytical studies of DNLS and its single impurity version 
12 1 . Concerning analytical results a little is known, only 
13 , [TJ, [lH . Recently it was proven for a rather general 
DNLS (even including disorder) that the energy spreads 
incompletely, provided that the norm which is a mea- 
sure of anharmonicity, is large enough [l5| . This proof is 
based on the existence of two conserved quantities, the 
energy and the norm. 

The main questions which now arise are: Does energy 
localization need the existence of more than one conser- 
vation law? Is there a sharp transition between com- 
plete and incomplete energy spreading? If so, what are 
properties characterizing such a transition? To explore 
these questions is our main motivation. Let us consider 
e.g. the /3-FPU model. In case that the energy delo- 
calizes completely one can linearize the equation of mo- 
tion at large times leading to the harmonic chain, which 
is exactly solvable. If it does not, it may converge to 
a DB. Since the amplitude of DB decays exponentially, 
one may linearize again, however outside the center of 
the DB. Idealizing this situation we consider a reduced 
FPU-chain with one anharmonic bond, only. The corre- 
sponding classical Hamiltonian for N particles with mass 
m and open boundary condition is as follows: 

N 1 c N ~ x 



2m 

n=l n=l,(n^M) 

+V(x M+1 - x M ) 



(1) 



where C > is the elastic constant and a n the equilib- 
rium length of the n-th bond. V(q) is chosen such that 
V has a single minimum at g m i n = with V"(q m in) = C . 
This rules out localized phonon modes of the linearized 
equation of motion. In addition we will assume that 
V"(q) is increasing with increasing \q — q m in\- H is trans- 
lationally invariant. After separation of c.o.m and in- 
troducing relative coordinates q n = x n +i — x n — a n and 
their conjugate momenta ir n , the harmonic part of Eq. ([T]) 
can be transformed to normal coordinates {Q^ , } and 
{Q% p *} for the left (1 < n < M - 1) and right 
(M + 1 < n < N) harmonic part of the chain. Skip- 
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ping the c.o.m. -energy, this yields: 



H = H\ 



harm 



(2) 



with -ff an h = — 7I"m + V(<1m) the Hamiltonain for the 
isolated anharmonic bond, the harmonic Hamiltonian 
Hharmi and Hi n t containing the interaction of the an- 
harmonic bond with the harmonic degrees of freedom 
(d.o.f.). Hamiltonian |J2J) is the only conserved quantity 
after separation of c.o.m. Since the equation of motion is 
linear in the harmonic d.o.f. these can be exactly elim- 
inated. This leads for N — >■ oo and M = 0{N) to the 
nonlinear integro-differential equation: 



1 



Q(r) + —V(q(r))- 



dr'k(r-r')^V'(q(T'))=0 (3) 



where the index M has been dropped for convenience. 
QUO) ee 0, P„ L (0) = and Q«(0) ee 0, P*(0) ee were 
chosen as initial conditions, r = u>ot is a dimensionless 
time and ujq = 2(C/m) 1 ' 2 the upper phonon band edge. 
The lower edge is at zero, due to translation invariance. 
The memory kernel is given by fc(r) = —ki(r) where 
fci(r) = J\(t)/t with J n the Bessel function of order n. 
Having determined for given initial conditions q(0) and 
q(0) a solution q(r) of Eq. ([3]) one obtains the harmonic 
nearest neighbor bond coordinates q n (r) from 



?»(t) = J dT'G lM _ nl (T-T')^V'(q(r')) , n^M 
o 

(4) 

with the Green function G n (r) = 2nJ2 n (r)/r. As initially 
localized excitation we choose q(0) = A and q(0) = 0. 
Use of a "velocity excitation" q(0) = 0, g(0) = B will 
not change our results qualitatively. With this initial 
condition in mind the conservation of the total energy 
implies that \q(r)\ < A for all r > 0. 

As well-known, elimination of a macroscopic number 
of d.o.f. induces dissipation. The frequency dependent 
damping constant 7(0;) follows from: 



j(u>) = lim 



1 -3 (k(u 
uj \ 



is) = 



Vi -uj 2 , M < 1 
, h>i 

(5) 

with uj measured in units of luq and k the Laplace trans- 
form of k(r). This exact result is obvious. For \u>\ < 1, 
i.e. for frequencies within the phonon band, the corre- 
sponding modes will be damped and consequently decay 
to zero, whereas all modes with frequency above that 
band will be undamped. If the anharmonic bond is iso- 
lated, i.e. the integral term in Eq. ([3]) is absent, q(r) 
will perform periodic oscillations with frequency Qq(A), 
depending on the amplitude A. Due to V"(q m i n ) = C 
it follows for ^4^0 that Oq(A) — > l/y/2 in units of uj . 
This frequency is within the phonon band. Since we have 
chosen a "hard" potential, i.e. dfip (A)/dA > 0, there will 



be a critical amplitude A^ such that Qq (A) touches the 
upper phonon band edge: 

O (4°)) = l . (6) 

Accordingly, one may speculate that for A < A^ the 
initial excitation will completely delocalize and will con- 
verge to a breather for A > A^. In the following we 
will chose a symmetric potential V{x)/C = \x 2 + jx 4 
for simplicity, x can be scaled such that the prefactor of 
the quartic term equals 1/4. In that case it is 

n (A) = ^V2 + AyK(-A 2 /(2 + A 2 )) (7) 

with K(m) the complete elliptic integral of first kind. 
Then Eq. © yields: 



Ai 0) 1.16715 



(8) 



In order to check the validity of our speculation above, 
we determine first the so-called limiting equation for the 
asymptotic solution (t) = liniA-»oo q(t + A.) [161 ]. The 
Laplace transform of Eq. ([3|) taking into account the ini- 
tial conditions can be solved for the Laplace transform 
q(z) of q(r) as function of q 3 (z). Transforming back to 
time regime yields: 



r 

q(r) = AJ (t) - J dr'J^r - r')<zV) 



(9) 



which is equivalent to Eq. ([3]) , as can be proven. For the 
pure harmonic chain, i.e. neglecting the nonlinear term, 
we obtain directly qhnrm{i~) = AJq(t), as is well-known. 
It is straightforward to derive the limiting equation: 



900(7-) 



T 



(10) 



Since q^ (r) is an asymptotic solution not possessing a 
relaxing component, its Fourier transform (foe (a;) can not 

have an absolutely continuous part c^(u>). If it would, 
its contribution q& (r) to q^ (t) would relax to zero for 
t — > 00 . Excluding a singular continuous component 
(which may occur for disordered systems at the mobil- 
ity edge), qoo(u>) must have a discrete support, i. c. 
goo(r) is either constant, periodic or quasiperiodic. If it 
is quasiperiodic, then there are at least two incommen- 
surate frequencies u)\ and u>2- The anharmonicity gener- 
ates Fourier modes with frequencies m\uj\ +miU)2- There 
exists an infinite number of integer pairs (mi, m^) such 
that m\W\ + rri2UJ2 is within the phonon band. There- 
fore, these modes are damped (cf. Eq. (J5J) and con- 
verge to zero. Accordingly consistent with our numer- 
ical results below, Eq. (fT0|) has two kind of solutions, 
only: A static one q^ atlc (T) = and a periodic one 
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r) with 2tt/tq > 1 in order 
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FIG. 1: n-dependence of A " from Eq. (TT3J for a — s,c. The 
dashed line represents ^ 1.16715 (cf. Eq. ©). 



to avoid an overlap with the phonon frequencies \uj\ < 1. 
Substituting g^ atic (r) = into Eq. HU]) yields the 
single solution q^ atlc (T) = 0. 

So far we have argued that two types of asymptotic 
solutions exist, a static and a periodic one. In order to 
investigate the existence of a critical amplitude A c we 
solve Eq. © iteratively. With the asymptotic behavior 
of J i we arrive at: 



q(r) ~ A^-[(r/r s y 

_i 

- (t/Tc) 2 CO: 
with relaxation times: 



sin | r — — 



r a (A) = 



]T(-iM Q) ^ 2 



a = s,c 



(11) 



(12) 



(3^ are given by n-fold integrals over products of J\ 
and J . Eq. (fTT|) with T a (A) from Eq. (fT2"| is a formal 
result for q(r) represented by a power series in A. It is a 
physical solution, only if the infinite sums in Eq. (j 1 2|) do 
exist. The critical value A c is such that this is true for 
A < A c . Then it is: 

A c = min{A c c \A c s n , = lim A^ , 

n — >oo 

1/2 



(13) 



An analytical calculation of these integrals seems impos- 
sible. Therefore it is done numerically which leads to A^ 
shown in Figure [T] up to n = 10. For n > 10 the numeri- 
cal errors become significant. This result gives evidence 
that A c is close to A c ^ . For A < A c the asymptotic time 
dependence of q(r) is similar to that of the harmonic 
solution AJq(t), however with a different phase and a 
renormalized relaxation time r re i(^4) = \J rf (A) + (A) , 
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FIG. 2: Top panel: A-dependence of q env (n) for n « 10 3 
(circles), 10 4 (plus signs) and 10 5 (crosses) time units. The in- 
set demonstrates the asymptotic behavior q env '(r») ~ A (solid 
line). Bottom panel: DB frequency Q, nwm (A) at n « 10 5 time 
units as function of A. The arrow indicates the critical value 
A { ° } from Eq. © and the dashed line A" um 1.181. The 
inset shows the asymptotic A-dependence Sl(A) ~ A (solid 
line). 



which diverges at A c . This behavior of T re i(A) follows 
from the divergence of the alternating sums (cf. Eq. §FI§ ) 
due to the quantitative difference of f3^ for n even and 
n odd, which also leads to the "oscillations" of An in 
Figure [T] According to Eq. (fTTj) , g(r) decays by an in- 
verse square root law, as also observed for the original 
/3-FPU chain [nj. 

In order to check these results and to access A > A c , 
we have solved Eq. ^ numerically up to r max = 10 5 us- 
ing an integration step of h = 0.05. Figure [2] depicts 
q env (n;A) for n « 10 3 , 10 4 and 10 5 where q env (T;A) 
is the envelope function of \q(r)\ for given A. With in- 
creasing Ti a clear sharpening of the transition is found 
at j4" um = 1.181, like for a second order phase transi- 
tion with finite size effects. A™ m differs from A { c 0) by 
about 1.2%. The frequency tt num (A) close to r max is 
shown in Figured] For A < A^ um we have Q num (A) ^ 1 
and for A > A" um it is well approximated by flo(A) 
for the isolated bond. However, for A above but close 
to A^ um the discrepancies are about 2%, whereas for 
A 3> A r c LUm they disappear. Whether the small devia- 
tion of A^ um and fl nam from A { c 0) and £1 {A), respec- 
tively, is genuine or stems from numerical inaccuracy is 
unclear. Hence it is not obvious that A c = A c ° . For 
A > yl" um the initial excitation indeed converges to a 
DB with frequency £l num (A). Figure [3] shows the numer- 
ically determined relaxation time T ie \(A) for A < A" um , 
and for A > A' c lurn the inverse modulation frequency 
2tt / uj moc \(A) of a modulation of the DB, which is ob- 
served numerically. For r — > oo the modulation am- 
plitude decays to zero. r re i has been determined from 
the criterion q env (T re i) = ^4/10. Both, r re i and 2-K/uo mo d 
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FIG. 3: Renormalized relaxation time r re i (circles) and modu- 
lation period 2ir/u> mo d (crosses) as function of A. The dashed 
line indicates A™ um = 1.181. The solid lines represent power 
law fits of T re ; and 2n/u>mod with exponents 0.61 and 0.87, 
respectively, which are supported by the log-log plots of the 
inset. 

seem to diverge at A" um by a power law with an expo- 
nent « 0.61 and « 0.87, respectively (see inset of Figure 
[3]). A power law divergence T re i(A) ~ (A c — A)~ 7 im- 
plies that /3„ ~ (A c ) _ "n~' 1_7 / 2 ^ for n — > oo. Whereas 
the exponential factor is strongly supported by our cal- 
culations the validity of the power law part can not be 
checked due to the limitation n < 10. 

Finally, we have analytically determined the moments 
ot) (r) = Eli(™-MM pot) (r), t= 1,2,3, . . . of the 

potential energy profile ei P ° t '' (r) in the thermodynamic 
limit. As a result we find 

T T 

m (pot)^ _ J dr\ J dTlKt(T - ti,t - T 2 ) 


x^V'Wt^V'Wt*)) . (14) 



Let us restrict to I = 2. K%{x^y) can be calculated ana- 
lytically and expressed by J\(x ± y) and Ji(x ± y)- Tak- 
ing into account the asymptotic expansion of J n we find 

m 2 ( T ) ~ 1-2 f° r T ~~ * °°! f° r a ^ ^- Hence, the en- 
ergy transportation is ballistic. This is expected since 
the transportation is within the half infinite left and right 
harmonic part of the chain. Using the profile of the ki- 
netic energy will not change these results. 

To summarize, based on combined analytical and nu- 
merical calculations of a reduced /3-FPU chain where the 
anharmonicity is restricted to a single bond we have pre- 
sented clear evidence for the existence of a critical ampli- 
tude A c which separates delocalization from localization. 
This demonstrates that a single conservation law is suf- 
ficient for such a transition. A™ um differs slightly from 

A^ . Therefore it is not clear whether A c coincides with 
A^ or not. Of course, no compelling arguments exist 
for their equality. In addition, the divergence of the iter- 
ation series is the mathematical origin of the transition 
at A c and leads for A < A c to a renormalized (due to an- 
harmonicity) relaxation time r re i which diverges at A c . 
The numerical solution suggests a power law divergence 
with exponent smaller then one. Above A c it yields the 
convergence towards a DB with frequency very close to 
£lo(A) of the isolated bond. Finally from the large r 
behavior of the second moment TO2(t) we find ballistic 
energy transportation for A < A™ um and A > A™ um . 
This proves that a divergence of W2(t) is not necessar- 
ily an indication of complete energy spreading, as it has 
been assumed for DNLS [ljj, supporting the conclusion 
in 0. 
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